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A characteristic of D-brane inflation is that fluctuations in the inflaton field can propagate at a 
speed significantly less than the speed of light. This yields observable effects that are distinct from 
those of single-field slow roll inflation, such as a modification of the inflationary consistency relation 
and a potentially large level of non-Gaussianities. We present a numerical algorithm that extends 
the inflationary flow formalism to models with general speed of sound. For an ensemble of D-brane 
inflation models parameterized by the Hubble parameter and the speed of sound as polynomial 
functions of the inflaton field, we give qualitative predictions for the key inflationary observables. 
We discuss various consistency relations for D-brane inflation, and compare the qualitative shapes 
of the warp factors we derive from the numerical models with analytical warp factors considered 
in the literature. Finally, we derive and apply a generalized microphysical bound on the inflaton 
field variation during brane inflation. While a large number of models are consistent with current 
cosmological constraints, almost all of these models violate the compact ificat ion constraint on the 
field range in four-dimensional Planck units. If the field range bound is to hold, then models with 
a detectable level of non-Gaussianity predict a blue scalar spectral index, and a tensor component 
that is far below the detection limit of any future experiment. 
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Understanding the physics of inflation [1] is one 
of the main challenges of fundamental physics and 
modern cosmology. Since string theory remains the 
most promising candidate for a UV completion of the 
Standard Model that unifies gauge and gravitational 
interactions in a consistent quantum theory, it seems 
natural to search within string theory for an explicit 
realization of the inflationary scenario. This search has 
so far revealed two distinct classes of inflationary models 
which identify the inflaton field with either open string 
modes (e.g. brane inflation [2-4], DBI inflation [5]), or 
closed string modes (e.g. Kahler moduli inflation [6], 
racetrack inflation [7], N-flation [8]). The brane inflation 
scenario [2, 3] in particular has received considerable 
theoretical and phenomenological interest. At the same 
time, precise cosmological observations [9] have made 
the detailed predictions of inflation testable. 

In this paper we study the general phenomenology 
of D-brane inflation models with arbitrary speed of 
sound. The theoretical ansatz is motivated by brane- 
anti-brane dynamics in warped spaces as described 
by the Dirac-Born-Infeld (DBI) action. Relativistic 
dynamics of the brane motion leads to deviations of 
the propagation speed of inflaton fluctuations from the 
speed of light. These effects generically result in a 
large level of non-Gaussianity of primordial fluctuations 
[10] and modify the inflationary single-field consistency 
relation [11]. We study these exciting observational 
signatures using a generalization of the inflationary flow 
formalism. This allows us to explore an ensemble of 
D-brane inflation models parameterized by the evolution 
of the Hubble parameter and the speed of sound as 
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polynomial functions of the inflaton field. 

The outline of the paper is as follows: In §2 we re- 
view the DBI mechanism of D-brane inflation and discuss 
its basic cosmological predictions. In §3 we discuss mi- 
croscopic consistency constraints that these models have 
to satisfy. In particular, we review and generalize the 
field range bound of [12]. In §4 we apply the inflation- 
ary flow formalism to D-brane inflation. We derive the 
generalized flow equations and use them to study the 
phenomenology of an ensemble of brane inflation models 
with arbitrary speed of sound. §5 contains a summary 
of our main results. We conclude in §6. An Appendix 
gives further details of the Monte- Carlo algorithm. We 
use natural units throughout the paper, where c = h = 1 
and Mp 2 = 8ttG. Interesting related work has appeared 
in [13] and [14]. 



2. DBI INFLATION 

2.1. Review of Warped D-brane Inflation 

Since a quantum theory of (super) strings is naturally 
defined in ten spacetime dimensions, any realistic de- 
scription of cosmology and particle phenomenology re- 
quires compactification of six of the space dimensions. 
Different compactification geometries lead to different 
low energy effective theories in four dimensions. In ad- 
dition, the extra dimensions have to be stabilized, since 
even a small time-dependence in the size and shape of the 
extra dimensions can manifest itself in a variation of fun- 
damental coupling constants and induce observable fifth 
forces. 

Recently, Kachru, Kallosh, Linde and Trivedi (KKLT) 
[15] provided a framework for stabilizing the extra di- 
mensions of type IIB string theory in the presence of 
background fluxes and non-perturbative effects (see [16] 
for a review of flux compactification). The flux back- 
ground fixes the shape (complex structure moduli) of 
the extra dimensions, but leaves the overall size (Kahler 
moduli) unfixed [17]. As [15] showed, the size of the 
compact space may be stabilized by the inclusion of non- 
perturbative effects, e.g. gaugino condensation on D7- 
branes. In addition to stabilizing the shape of the extra 
dimensions, the background fluxes lead to strong warp- 
ing of the spacetime. This warping provides an elegant 
mechanism to produce exponential hierarchies (a realiza- 
tion of the Randall- Sundrum scenario in string theory 
[18]) and has important consequences for the dynamics 
of brane motion in the warped space. 

The line element of a warped flux compactification of 
type IIB string theory to four dimensions takes the fol- 
lowing form 

ds 2 = h- 1 l 2 {y)g lxv dx^dx v + h 1 / 2 (y)g lj dy i dy^ , (1) 

where h is the warp factor and /i, v = 0, . . . , 3 and z, j = 
4, . . . , 9. Typically, the internal space will have one or 



more conical throats sourced by the background fluxes, 
i.e. regions in which the metric is locally of the form 



g ij dy i dy j = d X 2 ^X 2 ds 2 X5 



(2) 



for some five-manifold X5 which forms the base of 
the cone. If the background contains suitable fluxes, 
the metric in the throat region can be highly warped, 
h- 1 < 1. 

The KKLMMT scenario [3] refers to brane inflation 
[2] in a warped throat of a KKLT flux compactification 
[15]. In particular, most models of brane inflation are 
concerned with the dynamics of a mobile D3-brane that 
fills four-dimensional spacetime and is point-like in the 
compact space. Here, the inflaton field ip is identified 
with the geometrical separation between a D3-brane and 
an anti-D3-brane. The anti-brane is fixed at the tip of the 
throat (x = XiR ~ 0)? while the brane moves from large 
radius (x < Xuv) to small radius (xir)- 1 The dynamics 
of a D3-brane in the warped background (1) is governed 
by the Dirac-Born-Infeld (DBI) action 



M P 1l + 2V(X, if) 



(3) 



where 



p(x,<p) = -rV)vi - 2f{ V )x+r 1 {^)-v{ V ) . (4) 



Here ip = VX3X is the inflaton field, T3 is the tension 
of the D3-brane and / _1 (<£) = T^h~ l ((p) is the rescaled 
warp factor of the background spacetime. For a homo- 
geneous background X = —\g^ y d^(pd v (p « \(p 2 is the 
canonical kinetic term of the inflaton. The potential for 
the brane motion V((p) arises from moduli stabilization 
effects [3, 20] and the Coulombic brane-anti-brane inter- 
action. In the slow-roll limit fX <C 1, the DBI action 
(4) reduces to the following familiar form 



P(X,<p) = X-V(<p). 



(5) 



This slow-roll limit can be understood as the non- 
relativistic motion of the brane in the presence of a weak 
force from a flat potential. The relativistic limit of brane 
motion in a warped background is characterized by the 
parameter 7 (defined in analogy to the Lorentz factor of 
relativistic particle dynamics) 



1 



(6) 



Positivity of the argument of the square-root in (6) im- 
poses a local speed limit on the brane motion, Cp 2 < 
/ _1 ((/?) = Tsh -1 ^). The presence of strong warping, 



1 In this paper we consider the so-called 'UV-model', cf. [19] for an 
interesting alternative where the brane moves out of the throat. 
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Variable Description 


Notes 


7(0) 


Inverse sound speed 


equations (6) and (65) / Monte- Carlo 


H(</>) 


Hubble parameter 


Monte-Carlo 


e, 77, k 


Slow variation parameters equations (68) to (72) 




Scalar power spectrum 


We use fccMB — 0.02 Mpc -1 


Pt(k) 


Tensor power spectrum 


Pt = rP s 


r 


Tensor- to- scalar ratio 


equation (21); derived at /ccmb 


n s 


Scalar spectral index 


equation (19); derived at /ccmb 


n t 


Tensor spectral index 


equation (20); derived at /ccmb 


Inl 


NG parameter 


equation (27); DBI prediction. 


N e 


Number of e-folds 


dN e -Hdt 


X 


Radial coordinate 


equation (2) 


<P 


Canonical inflaton 


m 2 T 3 r 2 

T ■-> 


<\> 


Monte Carlo inflaton 


cf) — at the UV end of the throat. 


T 3 


D3-brane tension 


equation (35) 




Warp factor 


equation (1) and §3 




Warped brane tension 


f- 1 = T 3 /i _1 . Same as T(ip) in Ref. [13] 


9s 


String coupling 


9s < 1 


a' 


String scale 


a'Mp oc Vq- equation (36) 


M P 


4d Planck mass 


Mp — mpi/\J&K 


M, /v 


Flux on A- and 5-cycles 


integer quantum numbers 


N 


Five-form flux 


N = MK > 1 



NOTES. — A summary of variables and functions related to our description of DBI inflation and observables. 



h~ x <C 1, in the throat can make this maximal speed of 
the brane much smaller than the speed of light. When ip 
is close to this speed limit, then 7 is large. 

From the inflaton action (3) we find the homogeneous 
energy density in the field 



p = 2XV, X - V 

= (t-i)/- 1 ^, 

while the pressure is 

v = (i- 1 - 1 )r 1 -v. 



(7) 



(8) 



p and V source the dynamics of the homogeneous back- 
ground spacetime, ds 2 = — dt 2 + a(t) 2 dx 2 , as described 
by the Friedmann equations 



3Mj>H 2 = p 
2M 2 H = -(p + P). 



(9) 
(10) 



Accelerated expansion (a > 0) requires smallness of the 
variation of the Hubble parameter H = 9t lna, as quan- 
tified by the parameter 

e=-§t = l(l+w)<l, (11) 



where 



W= ~P = (7-l)/" 1 +V <-3- (12) 



From the expression for the equation of state parameter 
(12) we see that although the brane moves relativistically 
in the DBI limit, 7 ^> 1, inflation still requires that the 
potential energy V dominates over the kinetic energy of 
the brane (7 — ^T^h -1 . This is possible because the 
kinetic energy of the brane is suppressed by the large 
warping of the internal space, h~ x <C 1. 

As discussed in Ref. [21], the slow roll (non-relativistic 
brane motion) and the DBI (ultra-relativistic brane mo- 
tion) limits are connected continuously by an intermedi- 
ate regime where the relativistic effects are small, 7 = 
O(l), but non- negligible. 



2.2. Cosmological Observables 

To discuss the phenomenological predictions of the 
original DBI inflationary scenario [5] and its generaliza- 
tions [21], it is convenient to define the speed of sound 



dp 

dp 



P,x 



• ^A y,xx 7 



This is the speed at which fluctuations of the inflaton Sep 
propagate relative to the homogeneous background. In 
addition, we define slow variation parameters in analogy 
with the standard slow roll parameters for inflation with 
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canonical kinetic term 



XV. 



H 



H 2 M 2 P H 2 



e 

Cs 



(14) 

(15) 
(16) 



2.2.1. Perturbation Spectra 

To first order in the slow variation parameters, the ba- 
sic cosmological observables as a function of wavenumber 
k are [23] 



P.(k) 



Pt(k) 



1 H 2 



8tt 2 M|, c s e 



2 H 2 

n 2 M 2 



kc s —aH 



(17) 
(18) 



k=aH 



Scalar perturbations freeze when they exit the sound 
horizon kc s = aH, while tensor perturbations freeze 
when k = aH. Over a limited range of scales it is ap- 
propriate to parameterize deviations from perfect scale- 
invariance by the following spectral indices 



n s -l 



n t 



dP s 
dink 

dP t 
dink 



= — 2e — fj — k , 
= -2e. 



(19) 
(20) 



Notice the dependence of n s on the evolution of the speed 
of sound as captured by the parameter k. The tensor-to- 
scalar ratio in these models is 



r = — = 16c s e. 

-L. S 



(21) 



The dependence of (21) on the speed of sound implies a 
modified consistency relation 



> c s n t . 



(22) 



As discussed recently by Lidsey and Seery [11], equation 
(22) provides an interesting possibility of testing DBI in- 
flation. 

The standard slow roll predictions for the cosmological 
perturbation spectra are recovered in the limit c s — > 1, 
c, 0. 



2. 2. 2. Non-Gaussianity 

The non-trivial structure of the kinetic term in the ac- 
tion (4) leads to striking observational signatures of rel- 
ativists DBI inflation. In particular, it was observed [5] 
that this generically leads to a very large non-Gaussianity 
of the primordial fluctuations. 



Let us give a brief qualitative description of the phys- 
ical origin of this non-Gaussianity 2 before citing the re- 
sults of an exact computation [10]. Consider the unper- 
turbed kinetic term 



-'kin 



-r i v / w^=-^ 



and its first order variation under (p - 

£ (1) £kin = 1lf> 6(f>. 



7 

(p + Sep 



(23) 



(24) 



This indicates that self-couplings of cp are enhanced by 
factors of 7 arising from expansion of the square root 
of the DBI action. Non-Gaussianities come from the 
third-order interactions in £kin clue to the perturbation 
5ip(t,x). The leading effect can be estimated by con- 
sidering the ratio of the cubic perturbation to the mat- 
ter Lagrangian, S^jC^i n , to the quadratic perturbation, 
£kin, neglecting mixing with gravitational perturba- 
tions. Since matter self-interactions must dominate in or- 
der to obtain significant non- Gaussian features this pro- 
vides a rough estimate of the effect. In the large 7 limit 
the leading contribution to the non-Gaussianity parame- 
ter Jnl scales as [5] 



fNL OC (X 7 . 



5( 2 ) C 



kin 



(25) 



We observe that the magnitude of non-Gaussianities 
scales with the Lorentz factor 7. Observational con- 
straints on primordial non-Gaussianities therefore lead 
to interesting constraints on the magnitude of relativis- 
ts DBI effects. 

The non-Gaussianity of fluctuations in DBI inflation 
was estimated in the large 7 limit in [5] and computed 
more precisely for the general case in [10]. Observational 
tests of the non-Gaussianity of the primordial density 
perturbations are most sensitive to the three-point func- 
tion of the comoving curvature perturbations £. It is 
usually assumed that the three-point function has a form 
that would follow from the field redefinition 



c = c s 



(26) 



where Q g is Gaussian. The scalar parameter f^L then 
quantifies the amount of non-Gaussianity. It is a func- 
tion of three momenta which form a triangle in Fourier 
space. Here we cite results for the limit of an equilateral 
triangle. The general shape of non-Gaussianities in DBI 
inflation may be found in [10]. Slow roll models predict 
[25] fNL <C 1, which is far below the detection limit of 
present and future observations. For generalized infla- 
tion models represented by the action (3) with general 
pressure function V, one finds [10] 



Jnl 



35 
108' 



The following discussion parallels the treatment in Ref. [24] 
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where 



A = 



x 2 r,xx + lx 3 v tXXX 

XV tX + 2X 2 P ,xx 



(28) 



For the specific case of DBI inflation (4) the second term 
in (27) is identically zero [10] and the prediction for the 
level of non-Gaussianity is 



Inl = 



35 

108 \c 



:(7 2 -l). 



(29) 



Measurements of primordial non-Gaussianity therefore 
constrain the speed of sound during the time when CMB 
scales exit the horizon during inflation. This implies an 
upper bound on 7 from the observed upper limit on the 
non-Gaussianity of the primordial perturbations. Using 
the WMAP limit [26], -256 < f NL < 332 (95% confi- 
dence level), one finds 7 < 32. Ref. [27] forecasts con- 
straints on primordial non-Gaussianity from Planck of 
/atl 11 < 66.9 at the la level, which translates to 7 < 14. 

Using (29) the modified consistency relation (22) can 
be expressed completely in terms of observables [11] 



1 + 108 f 



(30) 



Although not inconceivable, it will be challenging to ex- 
perimentally test this unique prediction of brane infla- 
tion. For the standard slow-roll case with Sn t = — r, the 
combination of next-generation CMB polarization mea- 
surements of r and estimates of n t based on the direct de- 
tection of gravitational waves fails to test this consistency 
relation to better than 50% for even the most optimistic 
experimental scenarios [28]. 



3. MICROSCOPIC CONSTRAINTS 

The philosophy of this paper is to take maximal the- 
oretical guidance from the string theoretic origin of 
D-brane inflation models, but use a phenomenological 
Ansatz to capture the largest set of possible scenarios in 
a single framework. In particular, we allow considerable 
freedom on the functional form of the inflaton poten- 
tial V((p) and the background warp factor h(ip). Ulti- 
mately, both V and h should of course be derived from 
an explicit string compactification (see [22]). Here, we 
take the approach of studying a general set of functional 
forms for these quantities, restricted only by cosmolog- 
ical constraints (§2.2) and a minimal set of microscopic 
consistency requirements. In this section, we describe the 
microscopic constraints that we impose on our models. 



3.1. Warped Background 

3.1.1. Warped Geometry 

Schematically, the warp factor in equation (1) is deter- 
mined by the solution of the Laplace equation 



-V 2 h~\G 3 \ 



(31) 



where V 2 is the six-dimensional Laplacian and G3 pa- 
rameterizes 3-form flux which sources the warping. In- 
tuitively, equation (31) is just like the equation for an 
electrostatic potential on the compact Calabi-Yau space 
in the presence of a 'charge density'. The functional form 
of the warp factor h hence depends on how the 'charge 
density' \Gs\ 2 is distributed. If it is localized in one re- 
gion, then h will tend to decrease monotonically away 
from it. We call this the one-throat situation, h should 
reach some finite value (i.e. it should not diverge) at the 
bottom of the throat, since otherwise the warped string 
and brane tensions which scale as h~ l would vanish. In 
addition, to avoid a naked singularity h should not van- 
ish. 

The Klebanov-Strassler (KS) geometry [29] is an ex- 
plicit non-compact ten-dimensional solution to type IIB 
supergravity in the presence of background fluxes. The 
KS spacetime decomposes into the form (1) with the in- 
ternal space (2) given by a cone over X 5 = T 1 ' 1 . Far 
from the tip of the throat the KS solution is well approx- 
imated by AdS§ x T 1 ' 1 with the warp factor determined 
by the Green's function of (31) 



hAdS = 



where 



R 4 
("') 2 



TT 



Vol(X 5 ) 



X > Xm , 



N = MK. 



(32) 



(33) 



Here, Vol(X 5 ) parameterizes the dimensionless volume of 
X5 with unit radius. Typically, Vo\(X$) = (9(7r 3 ), e.g. 
Vo^T 1 ' 1 ) = ±ff^. g s < 1 and l s = Vc7 are the string 
coupling and the string length, respectively and xm is 
the minimal radial coordinate at the tip of the throat, 



HR/xm) 



2ttK 
3g s M 



(34) 



The integers K and M denote flux quanta on the A and 
B cycles of the throat. The exact KS warp factor is non- 
singular at the tip [29]. The AdS warp factor (32) forms 
the basis for many theoretical studies of brane motion in 
warped throat regions. However, other warped solutions 
are possible (and some are known), so in the spirit of our 
phenomenological approach and to retain maximal gen- 
erality we allow h(x) to be a free function subject only 
to minimal theoretical constraints. Of course, we appre- 
ciate that it is therefore not guaranteed that all (or even 
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most) warp factors we consider in this study have ex- 
plicit microscopic realizations. More detailed theoretical 
constraints on the functional form of h(x) than the ones 
presented in this section are beyond the scope of this pa- 
per. For a more complete theoretical understanding of 
the phenomenology of DBI inflation, such a theoretical 
study is essential. 

If one considers a scenario with two (or more) throats, 
corresponding to two localized 'charge' distributions in 
(31), then h will not be monotonic overall (similarly to 
the electric potential of two positive charges). It will 
reach a minimum between the throats ('charges') where 
h « 1. However, within each throat it should behave 
monotonically, as in the explicitly known examples of 
gauge/gravity duality like the KS solution. Let us remark 
that although we will focus our attention to single throat 
scenarios with monotonic warp factors, the methodology 
of the present paper is easily generalized to studies of 
multi-throat scenarios. 



3.1.2. Warped Brane Tension 

The dynamics of a D3-brane in a warped throat re- 
gion is determined by the brane potential V((p) and the 
warped brane tension / _1 ((/?) = T^h -1 ^). In the fol- 
lowing we relate the scale of the warped tension to mi- 
croscopic parameters. The D3-brane tension in four- 
dimensional Planck units is 



1 1 



(MA* 



M P (2tt) 3 g s \Mp) ' 



(35) 



where M~ 2 = a' defines the string scale. The string scale 
is related to the four-dimensional Planck mass (or the 
Newton constant G) via the (warped) compactification 
volume Vq u = J d 6 y^fgh 



MpY 



1 V? 



(36) 



Notice that the four-dimensional Planck mass increases 
as the (warped) volume of the internal manifold is in- 
creased. The warp factor in the throat region IrT 1 lies in 
the following range [h^, ^uvl> wnere 



V 1 



-87rK/3g s M 



(37) 



is the warp factor at the tip of the throat and h^y « 1 
defines the region where the throat is glued into a bulk 
space. Small string coupling and large K/M allow ex- 
ponentially large warping, i.e. small h^. Using / _1 = 
T 3 ft-\ we find 



(2tt) 3 5s 



M P 



I h 

(27T)" flf V if " 



(38) 
(39) 



To facilitate comparison with the analysis of Bean et al. 
[13], let us consider the same fiducial values for the string 
coupling and the string scale, g s = 0.1, (Mp/M s ) 2 — 
1000. This implies fixing the D3-brane tension to have 
the following value, T3 « 10 _8 Mp. Allowing the warp 
factor IrT 1 to be in the range h~ l G [10 -10 , 1] then implies 
fM% G [10 8 , 10 18 ] for the fiducial range of /. However, 
we note that by decreasing the string scale (increasing 
the volume of the internal space) and/or increasing the 
warping, one can make / much larger. In general, we 
treat the magnitude of / as a free parameter of the mod- 
els. 



3.2. Bound on the Field Range 

Since the inflaton field ip for D-brane inflation acquires 
a geometrical meaning, there exist geometrical restric- 
tions on the allowed range of cp. In particular, the finite 
size of the compact extra dimensions restricts the field 
variation of the canonical inflaton in four-dimensional 
Planck units. Naively, since cp oc x> one might imagine 
that the inflaton field range can be increased by simply 
scaling up the radial dimension of the throat. However, 
this increases the volume of the compact space, and by 
equation (36), changes the four-dimensional Planck mass. 
For the warped cones that form the basis of most ex- 
plicit theoretical models, the inflaton field range in four- 
dimensional Planck units in fact decreases as the radial 
dimension of the cone is increased. 

More specifically, recall from (35) and (36) that the 
four-dimensional Planck mass scales with the warped vol- 
ume of the compact space 



M 2 P = -(T 3 ) 2 V 6 w . 



(40) 



Using the conservative bound that the total warped vol- 
ume (which receives contributions from the throat and 
the bulk) is at least as big as the warped volume of the 
throat only, i.e. 



rXuv 

V 6 W > (^throat = Vol(X 5 ) / d X X 5 HX) , (41) 



finds 



For the cut-off AdS warp factor 



f A 

JAdS = —a • 

this implies that 



A = To,R A = - 



N 



2 Vol(Xs) ' 



Mp > j^v ■ 



(42) 



(43) 



(44) 
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and hence 



( <£uv\ 2 ± 



V M P J 



N 



(45) 



Since Acp < <^uv> these microscopic considerations imply 
the following bound on the total field variation during 
warped D-brane inflation [12] 



Mp) N ' 



(46) 



where for theoretical consistency the flux integer N = 
MK has to be much greater than unity. In all explic- 
itly known examples, super-Planckian field variations are 
therefore microscopically disallowed and brane-inflation 
models that predict Aip > Mp should be viewed with 
suspicion. 

For ultra-relativistic DBI inflation (/nl > 1) with a 
quadratic potential one can show that the normalization 
of the primordial scalar spectrum requires [12] 



N 



32 
3^ 



3Vol(X 5 ) 



> 10 8 Vol(X 5 ). 



(47) 



Since typically Vol(X 5 ) = 0(tt 3 ), observations therefore 
imply very large TV, allowing only very small field varia- 
tions. However, Ref. [12] also derived an upper limit on 
N in terms of the observational limits on non-Gaussianity 
and tensors 



27 

N<—rf NL <40. 



The limits (47) and (48) are clearly inconsistent unless 
Vol(X 5 ) is extremely small. 

The numerical analysis of Bean et al. [13] shows that 
the intuition gained from the explicit example of [12] ex- 
tends to the non-analytic, intermediate DBI regime. In 
particular, Bean et al. find that imposing the micro- 
scopic bound (46) dramatically reduces the parameter 
space of viable DBI models in the intermediate and ultra- 
relativistic regime. In this paper we study if this conclu- 
sion continues to hold when f(ip) is allowed to be a free 
function. 

For general f(ip), the total field variation during infla- 
tion is bounded by A(p < <^uv> or 



Mr) 



< 



■^throat 



where 



throat 



4Vol(X 5 ) 



* Vuv -W 



/ dip if 5 f 'fa). 



(49) 



(50) 



Typically, iVthroat > 1. In this paper we compute 7V t hroat 
by direct numerical integration of our output warp fac- 
tors f((p). Models that violate (49) cannot be embed- 
ded in a consistent string compactification. We find that 
most models discussed in this paper violate this condi- 
tion, which highlights the theoretical challenge of con- 
structing microscopically viable models of UV DBI infla- 
tion that are consistent with observations. 



3.3. Microscopic Bound on Tensors 

The field range bound has important implications for 
the expected level of gravitational waves from brane infla- 
tion models [12]. By the Lyth bound [30], the constraint 
on the evolution of the inflaton (46), or the generalized 
bound (49), is related to the maximal observable gravita- 
tional wave signal from inflation. Let us recall the argu- 
ment and generalize it to general speed of sound theories. 
Restricting to the homogeneous mode ip(t) we find from 
(14) that 



dip 
~Mp~ 



2e 



V. 



dN e 



and hence 



Acp 
~Mp~ 



pNcM 



■dN e 



x 



(51) 



(52) 



where Acp = <^cmb — ^end- Notice the non-trivial gen- 
eralization of the standard slow roll result through the 
factor V^x/j- For DBI inflation this factor happens to 
be unity, since V,x = 7, so that the Lyth bound remains 
the same as for slow roll inflation 



rcMB = WA^) ' 



(48) 

where 



pNcMB / 

/ diV e 

JN end \ 



eff 



fCMB 



1/2 



(53) 



(54) 



quantifies the support of the integral (52) and tqmb de- 
notes the tensor-to-scalar ratio evaluated on CMB scales. 
The value of 7V eff depends on the evolution of r(N e ) af- 
ter CMB scales have exited the horizon during inflation. 
In slow roll models this evolution is typically small (sec- 
ond order in slow roll), with observational data on CMB 
scales requiring iV e ff = (9(50 — 60), leading to the predic- 
tion of an unobservably small level of gravitational waves, 
rcMB <0.01 [12]. 



3.4. Implications for Relativistic DBI 

Lidsey and Huston [14] derived an interesting general- 
ization of the field range bound of [12] that is useful in 
the DBI limit (/jvl ^> !)• Here we briefly review their 
argument. Let Aip* be the field variation when observ- 
able scales are generated during inflation, corresponding 
to AN* < 4 e-foldings of inflationary expansion (this 
corresponds to the CMB multipole range 2 < £ < 100). 
Then the integral in the bound on the Planck mass (42) 
can be approximated as follows 



/ ' 



d W 5 /(^ > A^tfjA > (A^)V* • (55) 
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Here, we have bounded the integral by a small part of the 
Riemann sum, defined /* = /(</?+) and used A<y9+ < tp+. 
Equation (42) then becomes 



\M P 



< 



Vol(X 5 ) 



(56) 



Next, we note that the warped tension j~ x can be ex- 
pressed in terms of the scalar power spectrum P s , the 
tensor-to-scalar ratio r, and the non-Gaussianity param- 
eter fNL- f* can therefore be related to observables 



7T 

lCf 



Ptri 1 



1 



3/ 



NL 



and hence 



\Mp J 



< p* r 2 I 1 

16Vol(X 5 ) s * 



3/nl 



(57) 



(58) 



For slow roll models with ,/atl <C 1 this is not a very 
useful constraint. However, for relativistic DBI models 
with fjsfL ^> 1 the bound (58) is independent of /jvl- 
Since, P* ~ 2.4 x 10 -9 , r* < 0.5 (from observations) 
and Vo\(X$) = 0(tt 3 ) (from theory) we conclude that 
super-Planckian field variation is inconsistent with ob- 
servations. The Lyth bound (53) can now be written as 



Ay?* 
M P 



Substituting this into (58) we find [f 4] 



r* < 



32P* 



(AiV*) 6 Vol(X 5 ) 



(59) 



(60) 



The observed level of scalar fluctuations P* ~ 2.4 x 10 -9 
therefore implies that the tensor amplitude is unobserv- 
ably small for relativistic DBI models if the field range 
bound of [12] is applied (unless Vol(X5) is made unnat- 
urally small). We emphasize that the bound (60) does 
not apply to slow roll models since /nl > 1 has been 
assumed in its derivation. 

Interestingly, Lidsey and Huston also derived a lower 
limit on r for models of UV DBI inflation with /atl ^> 1 
and n s < 1 [14] 



> 



4(1 -n a ) 
V^7 



NL 



(61) 



The limits (60) and (61) are clearly inconsistent unless 
Vol(X 5 ) is very small, cf. (47) and (48). 

With the work of [12], [13] and [14] there is now a 
growing body of evidence that the best motivated the- 
oretical models of relativistic (UV) DBI inflation are in 
tension with the data if microscopic constraints are ap- 
plied consistently. In this paper we reach conclusions 
that are consistent with this and show that these prob- 
lems persist even if considerable freedom is allowed for 
the functional form of the brane potential V(ip) and the 
background warp factor f(ip). 



4. THE FLOW FORMALISM FOR D-BRANE 
INFLATION 

To study the ensemble of inflationary models speci- 
fied by the theoretical ansatz of the previous section, we 
adapt the inflationary flow formalism [31-33]. In this sec- 
tion we derive the generalized flow equations for brane 
inflation. In the next section we present our numerical 
results. 



4.1. Inflation in the Hamilton- Jacobi Approach 

Let us recall the Hamilton-Jacobi (HJ) approach to 
inflationary dynamics and apply it to warped brane in- 
flation. For a spatially flat FRW spacetime, the scale 
factor a(t) is determined by the Friedmann equation 

where p = (7 - l)/" 1 + V and 7 (0) = [1 - /(^)(/> 2 ]" 1/2 
follow from the DBI action (4). The inflaton 3 (f) obeys 
the following equation of motion 

where primes denotes derivatives with respect to <p. In 
the HJ-formalism the Hubble expansion rate H((j>) is con- 
sidered the fundamental quantity. The master equation 
relating <fr(t) and H((f>) follows from equations (62) and 
(63) 

*'« s f -f~**>5Bf- < M > 

From equation (64) one finds 

7(0) = ^l+Af{<t>)M P [H'(<t>)] 2 ■ (65) 



or 



-1 4[ff'] 



'12 



7 2 — 1 



Given H(<fr) and 7(0), the inflaton potential is 

[Hf 



V(4>) = 3MpH 2 - AM P 



7+1 



(66) 



(67) 



Notice the following important consequence of the 
Hamilton-Jacobi equations (66) and (67): For any speci- 
fied function H(cj)) and 7 (0), it produces a potential V((j)) 
and a warp factor f(<fr) which admits the given H((j)) and 
7(0) as an exact inflationary solution. 



3 In the following we use the variable <j> for the inflaton field rather 
than ip as in section 2, since we want to allow for the possibility 
that cj) = in our Monte Carlo simulation doesn't coincide with 
the tip of the throat. In fact, <f> = will be the UV end of the 
throat. The variables ip and (p are simply related by a linear 
transformation corresponding to this shift in origin. 
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4.2. Inflationary Flow Equations 

Analogous to the Hubble slow roll (HSR) parameters, 
we can write down a set of slow variation parameters for 
brane inflation that are defined in terms of derivatives of 
H and 7 with respect to <\>\ 



e(0) ee 

v(4>) = 
<4>) = 

e m = 

€ a(0) = 



2Mp 



7 (0) V H{4>) ) 
2Mp H"{<j>) 
7(0) W) ' 

7(0) H{4>)l{4>) 
f2M 2 P V (W_ 

U(0)/ v# 



(2M P V 



1 ^ +1 ff(0) 
H{4>) d^ 1 

1 ri w 7 (0) 
7(0) 



(68) 
(69) 
(70) 
(71) 
(72) 



where £ > 1. For notational convenience we define 77 = 
X A, £ = 2 A, p = 1 a and a = 2 a. Note that 77 is related to 
fj in (15) by fj = 2e — 2r] — k. The trajectories of these 
parameters are governed by a set of coupled first order 
differential equations. Using the relation 



2M 2 P H' 
7 ~W 



dN P 



-Hdt, 



we find 



de 

dNe 
drj 

dNe 
d £ X 

dN P . 



e(2e - 2ry + k) , 
r?(e + ^)+£, 
£k + £c-(£- 



A 



l)ry 



+ 



^+1 1 



and 



dn 

dNe 
dp 

dNe 
d £ a 

dK 



-k (2k ■ 
-2pK + 



e - rj) + ep , 



= — a 



{£ + l) K +(£-l)e 
l)ri 



+ e+1 a. 



(73) 

(74) 
(75) 
(76) 

(77) 
(78) 

(79) 



This set of differential equations defines the flow equa- 
tions for brane inflation. Notice that while the flow pa- 
rameters depend on 7, the flow equations do not depend 
on 7 explicitly. 

As pointed out by Liddle [34] , the flow equations have 
an analytic solution. Truncating the hierarchy of flow 
equations so that the last non-zero terms are M A and M a 
ensures that M+1 A=0 and M +1 a=0 at all times (along 
with all the higher order terms). From Eqs. (71) and 



(72), it then follows that higher derivatives vanish at all 
times: 



d0( M + 2 ) 
d0( M '+ 2 ) 



= o, 



0. 



(80) 
(81) 



We thus arrive at polynomials of order M + 1 and M' + 1 
in cj) respectively for the functions H((j>) and 7(0), 



H(4) = Ho 



1 + A ^W P 



+ A 



M+l 



\Mp 



M+l 



. (82) 



7(0) = 7o 



1 + Bl <Mp l+ - 



• • + Bm* 



M' + l 



Mr 



M' + l 



■ (83) 



So far, we haven't specified an initial value for <\>. If we 
truncate the series as above, any set of flow parameters 
spans an M + l dimensional space. However, the evo- 
lution equations define the flow parameters as functions 
of (j). Consequently, the set of distinct trajectories spans 
an M dimensional space, effectively fibering the space of 
initial conditions for the flow hierarchy. However, if the 
flow parameters are specified at <p = the ambiguity is 
removed. 

From the definitions of the flow parameters, the coef- 
ficients Aj and Bj can be written in terms of the initial 
values of the flow parameters 



and 



A x = 

A £+1 = 

B x = 

Bg+\ = 



\/ e o7o/2, 

(7o/2)' 
(£ + 1)! A[- x 



^o7o/ 2 
A, ' 

(7o/2)' 
(£ + 1)! A[- x 



Ao , 



a . 



(84) 
(85) 

(86) 
(87) 



The sign convention we choose is as follows. To define 
the direction of time with respect to the number of 
e-folds before the end of inflation 7V e , we choose N e to 
increase as one goes backward in time; i.e., dN e = —Hdt, 
so that dt > as dN e < 0. Furthermore, we choose y/e 
to have the same sign as H'(<j>). This is equivalent to 
choosing A\ > in our notation; the sign of A\ specifies 
in which direction the field is rolling. 

Note that in our present work f((j>) is not fixed to be a 
specific function. Instead, /(</>) is a derived quantity that 
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FIG. 1: The r — n s plane populated by numerical models of DBI inflation with M — 5 and M' — 1 in equations (82) and (83), 
combining simulations with general and small 70 priors (left panel) and by numerical models of standard slow-roll inflation 
from Kinney [36] (right panel). 



is determined through equation (66) via Monte-Carlo de- 
scriptions of 7(0) and H((j>) . Similarly, V((j>) is also a de- 
rived function determined through H((j)) and 7(0) using 
equation (67). In this sense, our work differs significantly 
from the recent work of Bean et al. [13] in that they 
fix both V(<fr) and f(<fr) to particular (well- motivated) 
forms. We allow a large range of warp-factors through a 
specification of a general sound speed captured by 7(0), 
and a large range of dynamics captured by H(<fr). Thus, 
the observable distributions shown in Bean et al. [13] 
for quantities such as the scalar spectral index and the 
tensor-to-scalar ratio are under a fixed theoretical model 
whose internal parameters (such as the normalization of 
the AdS warp factor) are allowed to vary. This can be 
considered a 'top-down' approach. Our work adopts the 
complementary 'bottom- up' viewpoint that, while some 
(or even many) of our warp factors may not have mi- 
croscopic realizations, it is still interesting to investigate 
the phenomenology of the set of warp factors, potentials, 
and observables that are generally allowed by existing 
cosmological data. 



4.3. Monte-Carlo Algorithm for DBI Inflation 

Now we shall outline the algorithm used to produce the 
numerical results of §5; technical details of the specific 
implementation are described in Appendix A. 

The algorithm is designed to capture the dynamics of 
DBI inflation from the time the brane enters the mouth 



of the throat and inflation begins, until inflation ends 
(either through a tachyonic transition near the tip of the 
throat or due to slow roll ending). Because we explicitly 
define the direction of motion as the brane moving into 
the throat (0 decreases), the algorithm as described be- 
low only applies to monotonic warp factors (single throat 
scenarios). However, it is straightforward to generalize 
the method to include the brane moving out of the throat 
((/> increases) and thereby model multi-throat scenarios 
with non- monotonic warp factors; we plan to explore this 
possibility in future work. 

A key feature of the algorithm is a well-defined physical 
"location" for (j) = 0, corresponding to the mouth of the 
throat, 0uv- Here and in the following, the subscript "0" 
denotes evaluation at = 0. 

There are three phenomenological classes of DBI in- 
flation that we would like to consider in this work. We 
classify models by their value of 7 at CMB scales: 

1. slow roll DBI: 1 < 7cmb < 1.1 

2. intermediate DBI: 1.1 < 7cmb < 10 

3. ultra-relativistic DBI: 7cmb > 10. 

We are particularly interested in the limit of the DBI 
classes above which are consistent with the microphys- 
ical bound (49) on the field range. In order to sample 
these different regimes efficiently, we consider different 
combinations of priors on 70 and eo in our Monte-Carlo 
simulations, as follows: 
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FIG. 2: Cosmological observables evaluated at /ccmb for a simulation applying the general 70 prior with M = 5 and M' — 1 
in equations (82) and (83). The three regimes of DBI inflation are shown with color-coded points: slow roll DBI (black), 
intermediate DBI (orange) and ultra- relativistic DBI (blue). The individual plots from top-left in clockwise direction are: (a) 
the scalar-to-tensor ratio r vs. scalar spectral index n s , (b) spectral index n s vs. running of the spectral index a s — dn s /d\nk, 
(c) r vs. Acj) = 0o — 0end, and (d) r vs. non-Gaussianity parameter /nl- In (a), we show the lower limit on r for ultra- relativistic 
DBI models with r > (1 - n s )/8 from [14] as a dashed line. In (d), the limit from WMAP, |/ NL | < 332 [26] is indicated by a 
vertical dashed line. In (c), the relation between r and Acj) can be described approximately as r oc (A0/M P ) 2 . This follows 
from equation (53); integration from iVo to iY en d instead of from Ncmb to iV en d just results in a different normalization iY e ff. 
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FIG. 3: Same as Fig. 2, but with a small 70 prior. Note that while in Fig. 2(a), the tensor/scalar ratio for ultra-relativistic 
DBI models satisfied a general lower bound, this is not the case for non-relativistic models. 



• Small 70 initial conitions: We assume that the initial 
speed of the brane is very small, so 

70 = 7 (0) = 1 + A, (88) 

for A < 1 and draw log 10 70 randomly from the 
narrow flat prior [0, 10 -5 ]. This prior is effective 
at sampling the slow roll and intermediate DBI 
regimes well. 

• General 70 initial conditions: We relax the above as- 



sumption on the initial speed of the brane, drawing 
log 10 7o randomly from the broad flat prior [0,1]. 
This prior is effective at sampling the intermediate 
and ultra-relativistic DBI regimes well. 

• General eo initial conditions: We allow a key observ- 
able quantity, the tensor-to-scalar ratio, to take a 
wide range of values by drawing from the broad flat 
prior log 10 eo G [—10,0]. Unless specified otherwise, 
this is the standard prior we will use in §5. 
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• Small eo initial conditions: As we will find below, the 
general eo prior generates DBI models that are con- 
sistent with the current cosmological constraints, 
but which strongly violate the microphysical bound 
(49). In order to quantify the limit on the tensor- 
to-scalar ratio at which the field range bound is 
satisfied for the different classes of models, we draw 
from the narrow fiat prior log 10 eo G [—21, —9]. 

The rest of the algorithm is identical for all the cases 
that we consider. We Monte-Carlo over Pq = P s (4> = 0), 
the initial amplitude of the power spectrum of scalar den- 
sity perturbations at (j) = and its value at the CMB 
scale P s (&cmb)- The initial values at the UV end for a 
limited number of the slow variation parameters in equa- 
tions (82, 83), which one expects may be constrained by 
current and future cosmological data, are also picked ran- 
domly. The priors for these initial conditions are specified 
in §A.l. The initial conditions at <p = fix fo and Hq, 



fo 



H n = 



rt- 1 m- 4 



1 



(89) 
(90) 



We numerically evolve the flow equations forward into 
the throat (</> decreases) until we find a match to the 
scalar power spectrum amplitude at CMB scales: 



7CMB ( HCMB \ ' 
87T 2 e C MB \ M P ) 



P S (kcMB) 



(91) 



The way that the matching condition at CMB scales is 
implemented is described in detail in §A.2. Once this con- 
dition is satisfied, we can compute the observable quanti- 
ties by linking the value of the infiaton field <p to the co- 
moving wavenumber k(fa) of cosmological perturbations, 
as follows. 

Without loss of generality, we pick some fiducial phys- 
ical scale kcMB to correspond to <\> = 0cmb- Then, with 
our sign convention, <\> > </>cmb corresponds to scales 
larger than &cmb , and (j) < (/>cmb corresponds to smaller 
scales. The wavenumber of cosmological perturbation 
modes is associated with a value of <\> through [21, 23] 



dNe 
dink 



1 



1 — e — k 



so 



d^ 

a M P 

d\nk 



1 1 



(92) 



(93) 



where the last expression follows from equation (73). 
Hence we associate a wavenumber with a value of (j) by 
solving Eq. (93). 

Instead of the amplitude of the power-spectra at each 
/c, the observables are widely described in terms of the 
power-law variables. To second order in slow variation 



parameters the spectral indices of the scalar and tensor 
perturbations are [13, 21] 



1 



n t 



dP s 
d\nk 

dP t 
d\nk 



(l + e + /t)(-4e + 2r ? -2 K ), (94) 

(») 



1 



We also consider the variation with scale, or 
of the spectral indices 



running 



a. 



a t 



dN e 



o dr l 



(i 



dn s 
d\nk 

dn t 

d\nk (1-e-ft) 2 



dN e 



O de 

dN e 



(96) 
(97) 



In our simulations we compute the cosmological ob- 
servables P s (k), Pt(k), and Jnl, as well as the various 
power-law variables, at &cmb = 0-02 Mpc -1 [38, 39] us- 
ing the k(fa) equation (93). Cosmological scalar modes 
freeze as they exit the sound horizon, c s iJ _1 = (jH) -1 , 
while tensor modes freeze when their scale exceeds the 
Hubble radius, H -1 . Therefore, /ccmb corresponds to 
different values of the field (j) s and fa for the scalar and 
tensor modes respectively, and this difference is taken 
into account when computing observables. For the ten- 
sor/scalar ratio, we calculate P t {fa)/ P s {4> s ), rather than 
using the analytic expression for r. 

In addition, throughout the evolution, we impose a 
set of constraints to enforce consistency of the physical 
picture, which we describe in detail in §A.2. 

Now we will consider the implementation of the end 
of inflation in our numerical models, leaving a detailed 
description until §A.3. In most models brane inflation 
ends via a tachyonic instability as the separation between 
a D3-brane and an anti-D3-brane becomes comparable 
to a string length. This correlates with the warp factor 
reaching a minimal value since the anti-brane minimizes 
its energy at the tip of the throat where the warp factor 
is smallest. Here, we consider two distinct scenarios: 

1. Inflation ends by tachyonic instability. 

2. Inflation ends before the tachyonic instability sets 
in, i.e. e — > 1 before the D3-brane comes within 
one string length of the anti-brane. The tachyonic 
mode then only serves to remove the inflationary 
energy density, so that the vacuum has zero energy 
after reheating. 

One more ingredient is needed to complete the descrip- 
tion of our algorithm. We found that simply Monte- 
Carlo'ing the priors for initial slow-roll parameters as out- 
lined in §A.l produced large numbers of models which, 
while satisfying all our physical constraints, had observ- 
ables which were very different from the current cos- 
mological data. This is because the evolution of flow- 
equations is very sensitive to the initial conditions due 
to the rapid variations in the function 7(0), which when 
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combined with H(<fr) leads to rapid variations in f (</)). 
While the goal of this present work is not to make a de- 
tailed comparison of the model with the data, we are nev- 
ertheless primarily interested in the properties of gener- 
alized DBI models whose evolution histories are broadly 
consistent with current observations. 

In order to preferentially increase the population of 
such models in our simulations, we apply the following 
selection mechanism: a Metropolis-Hastings algorithm 
(which is used in standard Markov Chain Monte Carlo 
techniques) is implemented with a "penalty function" 
taken from the scalar P(k) results of Fig. 10 (bottom 
right panel) from Ref. [39]. This figure contains a scalar 
P(k) reconstructed from WMAP 3 year data [41, 42] 
and the SDSS galaxy power spectrum [43], under the 
assumption that the primordial fluctuations are seeded 
by the standard single field slow roll inflation mecha- 
nism that additionally satisfies a minimal "sufficient e- 
folds" requirement that solves the cosmological flatness 
and horizon problems [44]. Since we expect, in general, 
that the present model, which contains more parame- 
ters, will be less well constrained by the current data, 
the Metropolis-Hastings algorithm is run at a high "tem- 
perature" . In practice, this means that the penalty func- 
tion, a least-square statistic, uses double the 95% CL 
error of the figure as a 1-cr Gaussian error, allowing the 
DBI model considerably more freedom to deviate from 
the mean than the single field slow roll case. This works 
extremely well in practice to find models broadly compat- 
ible with the data. It is straightforward to incorporate 
this module into standard parameter estimation codes to 
do a direct comparison with cosmological data, and we 
will present such results in future work. 



5. RESULTS 

In Figures 1 to 11, we present our results for observ- 
ables at CMB scales as well as the functions associated 
with the underlying models. We classify DBI models 
by their value of 7 at CMB scales: 1 < 7cmb < 1-1, 
1.1 < 7cmb < 10, 7cmb > 10 define the slow roll (black), 
intermediate DBI (orange) and ultra-relativistic DBI (blue) 
regimes, respectively. 

Different combinations of the priors on 70 and eo de- 
scribed in §4.3 have been applied to each figure. All fig- 
ures except Fig. 11 use the general eo prior; the latter has 
the small eo prior applied. As specified in the figure cap- 
tions, sometimes we will show the general and small 70 
prior simulations separately to emphasize different points 
of interest, and sometimes we will combine these simu- 
lations. While in general the plots with combined priors 
have different measures on models in each simulation, 
this does not matter for our purpose, which is to investi- 
gate the functional parameters and predictions of these 
models qualitatively. 

In addition, all figures except Figs. 9 and 10 show sim- 
ulations run with M = 5 and M' = 1 in equations (82) 



and (83); the latter used M = 2 and M' = 1. In the 
case with six H(cj)) parameters, we found that all of the 
models had inflation ending with tachyonic transitions. 
This is because models with significant higher order 
derivatives of H(<fr) which end inflation with e = 1 
have a strong tendency to provide insufficient e-folds of 
inflation to satisfy our requirements, and are therefore 
rejected from our simulations. In order to get a fraction 
of models where inflation ended with the end of slow roll, 
e = 1, we had to reduce the number of H((j>) parameters 
to three; even in this case, tachyonic transitions were the 
dominant mechanism for the end of inflation. In Figs. 9 
and 10 we show the properties of only the models where 
inflation ended with e = 1. 

Here we highlight and comment on what we consider 
to be the most important results of this study. Further 
details can be found in the figures and the figure captions. 

Fig. 1 compares the distributions of points in the n s 
vs r plane from a Kinney-style single-field slow roll sim- 
ulation [32] and a DBI simulation. Note that the algo- 
rithm used to create the slow roll simulation is signifi- 
cantly different from ours; among other dissimilarities, 
it makes use of both forward and backward integration 
of the single-field slow roll flow equations, while our al- 
gorithm only uses forward integration of the generalized 
flow equations. The slow roll simulation is shown here 
because it adopts a similar philosophy that "fundamen- 
tal" quantities of interest which are not strictly fixed by 
theory, such as the potential, should be Monte-Carlo'd 
rather than fixed to specific functions when investigating 
the range of observable properties predicted by a given 
model. More detailed inferences should not be made from 
this comparison because of the differences in algorithms. 

A paucity of points is apparent in the DBI simulation 
as \n s — 1| — > 0. We explain this as follows: For per- 
fect AdS warp factors f(<fr) oc (0 — (j) en d)~ A and a linear 
Hubble parameter H((j)) oc (cj) — (j) e nd) one may show an- 
alytically that n s = 1 to all orders in inflationary slow 
roll parameters [21]. However, it has also been shown 
that non-linear corrections to H((j)) and deviations of the 
warp factors f((j)) from the AdS limit lead to deviations 
from scale invariance [40]. This is what we are seeing 
here. Our warp factors are generally flatter than AdS for 
the regime of cosmological interest, and H"(<fr) is often 
non- negligible. Below we explain why we believe our al- 
gorithm preferentially selects non- AdS warp factors. The 
same argument then explains the lack of models leading 
to perfect scale invariance. 

Fig. 2 and 3 present the results for the basic cosmo- 
logical observables evaluated on CMB scales for general 
70 prior and small 70 prior, respectively. As illustrated in 
Fig. 2(a), which shows r vs. n s , models with detectably 
large values for the non-Gaussianity parameter Jnl sat- 
isfy a lower limit on the tensor-to-scalar ratio r if the 
scalar spectrum is red n s < 1. This bound was derived 
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FIG. 4: Slow-roll inflation (left) and DBI (right) consistency relations, combining the simulations with general and small 70 
priors shown in Figs. 2 and 3, respectively. The left panel shows r vs. n*. In standard slow-roll inflation r = — Snt (dashed 
line), but as shown, DBI inflationary models depart significantly from this relation except in the case of slow roll DBI models 

(black points). In the right panel, we show r [l + ^/nl} 1 ^ vs. — 8n t describing the DBI consistency relation; all models 
satisfy this relation. 




FIG. 5: A consistency check of the field range bound (49) for the general 70 simulation shown in Fig. 2 (left) and the small 70 
simulation shown in Fig. 3 (right), with A0 2 /Mj> vs. the expression on the right hand side of Eq. (49), assuming Vol(Xs) = tt 3 . 
The microscopic bound requires that A(/> 2 /Mj> be smaller than the right side of eq. (49), and all models in these simulations, 
which had a general eo prior applied, violate this bound. 
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FIG. 6: rfl/l vs. 1 — n s for general 70 (left, Fig. 2) and small 70 (right, Fig. 3) simulations of DBI models. The black line is 
the lower limit on tensors for relativistic DBI models [14]: r > 4 (j~ ns ^ . While this relation between n s and r is satisfied by 

ultra-relativistic DBI models, this relation is not satisfied by the intermediate and slow roll limit DBI regimes. 



in [11, 14] and may be written as 

4 



r > 



(98) 



Notice that the bound (98) only applies if /nl > 1, as 
can be seen from the r vs. n s plot in Fig. 3(a). Since 
current CMB and large-scale structure data are strongly 
indicating a red scalar spectrum [9], Fig. 2 then implies 
the exciting prediction that DBI models with detectably 
large non-Gaussianities correlate with detectably large 
tensors. Unfortunately, this conclusion ignores an im- 
portant theoretical consistency constraint that should be 
imposed on any viable model. As Ref. [12] showed, the 
inflation field range in warped D-brane inflation is geo- 
metrically limited, and in explicit models can be written 
in the following simple form 



Acf) 2 



(99) 



where N ^> 1. By the Lyth bound [30] this impossi- 
bility of super-Planckian field variation during inflation 
implies a strong upper limit on r. For the generalized 
warp factors studied in this paper the field range bound 
is written in the generalized form (49). Fig. 5 illustrates 
this bound for the models shown in Fig. 2 and 3. The 
plots show rather dramatically that all models in Fig. 2 
and 3 violate the microscopic bound of [12]. 

To find models that satisfy the field range bound (49) 
we ran separate simulations with small eo prior. The re- 



sults are illustrated in Fig. 11. Here, the y-axis rep- 
resents the compactification constraint on the warped 
throat volume (42) written in the following form 

This condition has to be imposed as a consistency con- 
straint on all models. Models that violate (100) corre- 
spond to a mismatch between the field range required by 
the specific inflationary model and the field range allowed 
by the compactification. We find that (100) can only be 
satisfied if the field range from IR (</> en d) to UV (0o) is 
very small (in Planck units). This is because the magni- 
tude of / is constrained to be large by the normalization 
of the scalar spectrum, so the size of the integral can't 
be made small by making / small. Instead the integral 
can only be small if the range of integration is small. By 
the Lyth bound this small range for the field variation 
corresponds to a very small tensor signal. 

From Fig. 11 we see that theoretically consistent inter- 
mediate and ultra-relativistic DBI models satisfy a very 
strong upper bound on tensors: r < 10 -15 . This limit 
is consistent with the bound (98) since the second plot 
in Fig. 11 shows that all consistent models with large 
/nl (so that (98) applies) have a blue spectrum n s > 1. 
These models could therefore be falsified in the near fu- 
ture if the current indications of a red spectrum gains in 
statistical significance to be more than 3<r. 

We emphasize, however, that the bound of r < 10 -15 
is not absolute and in particular does not apply to slow 
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FIG. 7: Input and output functions for the general 70 simulation from Fig. 2. From top- left in clockwise direction are: (a) 
Ps(k), where the dashed lines are the 68% and 95% CL observational constraints from Ref. [39], (b) /(</>), where the dashed 
line illustrates the shape of the AdS warp factor with f(<j>) oc ((f) — en d) _4 , (c) V ((/)), and (d) 7(0). 



roll models with small Jnl- The upper limit on r for 
consistent slow roll models (black points) in Fig. 11 de- 
pends on the lower limit of the prior on 70. For very 
small 7, larger r becomes consistent with the field range 
bound. (Similar conclusions were obtained by Bean et al 
[13]. In particular, they showed that only extreme slow 
models with 7cmb — 1 < 10 -7 survive the bound of [12]). 
However, the theoretical considerations of [12] still apply, 
which predict an ultimate upper limit on r for slow roll 
models. This limit is far below the detection sensitivity 
of any realistic future experiments. 

We now comment on the shapes of the empirical f(<f>) 
derived from our simulations as compared to the analytic 
AdS warp factor f(<j)) oc <p~ A studied in the literature in 
connection with DBI models. In our simulations with 6 
H((j)) parameters, many of the recovered f(<fr), as seen in 
Figs 7(b) and 8(b), exhibit power-law behavior similar to 
the AdS warp factor, but are somewhat shallower than 
the AdS case. Further, there are some models which start 
out with power-law behavior and flatten out towards the 
end of inflation. In the case of our simulations with 3 



H((f)) parameters, where we consider models where infla- 
tion ends with e = 1, the derived f(<p) shapes are ex- 
clusively of the latter kind. Our work shows that, qual- 
itatively, there are some differences between the shapes 
of warp factors found through the empirical process de- 
scribed in our algorithm and the analytical considerations 
so far. However, given the large parameter space explored 
in our simulations, there was no a priori expectation that 
we would recover this particular function. Therefore, it 
is encouraging that the derived functions are not, in fact, 
dramatically different from the analytic form. Given that 
our prior is broad enough to allow warp factors to be 
generated in a Monte-Carlo fashion, two possibilities are 
that either the shallower functional forms are more likely 
to be encountered in the prior, or that the penalty func- 
tion used in our Metropolis-Hastings algorithm prefers 
shallower functions over steeper ones. In a future study, 
we plan to make a full investigation of the likelihood of 
the analytic AdS warp factor compared to the derived 
empirical warp factors, in order to test whether the data 
in fact prefers shallower functions for f ((/)). 
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FIG. 8: Same as Fig. 7, but for the 
6. CONCLUSION 

In this paper we studied the basic phenomenology of D- 
brane inflation models with general speed of sound. We 
developed a general Monte-Carlo formalism for studying 
models with generalized warp factors and inflaton poten- 
tials. Most models show large deviations of the speed of 
sound from the speed of light resulting in exciting ob- 
servational signatures. Non-Gaussianities are typically 
large enough to be observable and the standard slow 
roll consistency relation is violated at a non-negligible 
level. However, we also showed that most of these phe- 
nomenological models cannot be embedded into a con- 
sistent string compactification. In particular, the vast 
majority of models violates the field range bound of [12]. 
UV models of DBI inflation which obey the field range 
bound are of two distinct types: 4 




□.001 0-Q10 0,100 1,000 1 0.00C 

E H(flfii-} pcK-n meters, Small y 
1 00F — 1 ~^ ! " " J "J 




0.001 0.010 0,100 1,000 1 0.00C 

small 70 simulation from Fig. 3. 

1. Slow roll models with unobservably small /nl- 

2. Relativistic DBI models with observable Jnl, but 
blue scalar spectrum, n s > 1. 

Let us therefore imagine that future observations yield a 
firm measurement of a red scalar spectrum, i.e. a sta- 
tistically significant detection of n s < 1 after marginal- 
izing over (or detecting) r and dn s /d\nk and the "late- 
time" cosmological parameters. What would this imply 
for brane inflation in general and the DBI limit in par- 
ticular? As indicated above, such an observation of the 
scalar spectrum would effectively rule out relativistic DBI 
models, without having to perform any measurements of 
non-Gaussianities. This illustrates that cosmological ob- 
servations are quickly becoming precise enough, so that 
many theoretical ideas are highly constrained. DBI in- 
flation is highly falsifiable and might in fact very soon 
be ruled out by the data. 5 This conclusion motivates 



4 In addition, there is the 'IR model' [19] which requires a separate 
analysis. 



5 To be precise, we emphasize that this statement at present only 
applies to realizations of the DBI mechanism using D3-branes 
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FIG. 9: Observable parameters at /ccmb for a general 70 simulation with M = 2 and M' — 1 in equations (82) and (83) 
respectively, showing models where inflation ends with e — > 1. The individual figures from top- left in clockwise direction 
are: (a) the scalar-to-tensor ratio r vs. scalar spectral index n s , (b) spectral index n s vs. running of the spectral index 
a s = dn s /d]nk, (c) r vs. n*, with a dashed line showing the standard slow roll expectation that r — —8n t , (d) A0 2 /M|> vs. the 
expression on the right hand side of Eq. (49) assuming Vol(Xs) = 7r 3 , again showing that these models violate the microscopic 
bound on the field range, (e) r vs. A</> = 0o — en d, and (f) r vs. non-Gaussianity parameter /nl, showing the limit from 
WMAP, |/nl| < 332 [26] with a vertical dashed line. The equivalent small 70 simulation gives qualitatively similar results. 



considering slow roll models of brane inflation. However, 
these models suffer from a version of the supergravity 
eta-problem [3]; moduli stabilization effects typically in- 
duces large corrections to the inflation mass and inflation 
cannot occur. Given the prospect that relativistic DBI 
inflation via D3-brane motion on Calabi-Yau cones might 
soon be ruled out by the data, it becomes important to 
understand whether the eta-problem for slow roll brane 
inflation can be overcome [22, 51, 52]. 



on Calabi-Yau cones. In this sense our conclusions do not repre- 
sent model-independent constraints on the DBI mechanism. In 
particular, all our statements are restricted to the concrete re- 
alizations of DBI inflation that we described in §2. Given that 
the strong constraints on these models are mostly driven by the 
geometrical field range bound of [12], it would be very interest- 
ing to find generalizations of this class of models that evade the 
bound. Concrete efforts to construct such scenarios are under 
way [49, 50]. 



We are living in the fortunate age where the cosmo- 
logical data are becoming precise enough to significantly 
constrain theoretical models of the early universe. 
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FIG. 10: Input and output functions for the models in Fig. 9. 
equivalent small 70 simulation gives qualitatively similar results 
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APPENDIX A: MONTE-CARLO ALGORITHM 
FOR DBI INFLATION 



In §4.3, we gave a broad outline of the Monte-Carlo 
algorithm that we use to produce the numerical results 
of §5. In this Appendix, we describe its implementation 
in detail. 



The sign convention is tied as follows to the physical 
picture of a brane evolving from the UV mouth of the 
warped throat towards its IR tip. As time increases, the 
brane falls into the throat towards the end of inflation. 
We will call this direction forward. d<p < as the brane 
moves forward. N e = at the end of inflation by conven- 
tion, and large at (j) = 0; i.e. 7V > Nqmb > ^V"end where 
-Wend = 0. In practice, we set Nq to some large value, and 
once we find N en( ^ by evolving forward, we just translate 
the zero point so that N en( ± = 0. 



1. Initial Conditions 

First, since we define (j) = be the UV end of the 
throat, then by definition, 

h£, = hr\ti) = \, (Al) 

and 

/o" 1 = /" 1 (0)=l3. (A2) 

As stated in the main text, we apply several different 
types of priors on 70 and eo, as follows: 

• General 70 initial conditions: We draw log 10 70 ran- 
domly from the broad fiat prior [0, 1]. 

• Small 70 initial conditions: We draw log 10 70 ran- 
domly from the narrow fiat prior [0, 10 -5 ]. 

• General eo initial conditions: We draw log 10 eo ran- 
domly from the broad fiat prior [—10,0]. 

• Small eo initial conditions: We draw log 10 eo ran- 
domly from the narrow fiat prior [—21,-9]. 
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FIG. 11: (a) Tensor /scalar ratio r (left) and (b) scalar spectal index n s (right) vs. the microphysical constraint on the left hand 
side of Eq. (100) assuming Vol(Xs) = 7r 3 , combining general and small 70 priors with M — 5 and M' — 1 in equations (82) 
and (83). Unlike previous simulations shown, these models have a small eo prior applied. The points below the horizontal 
line are DBI inflationary models allowed by the microphysical bound while models above this line violate the bound and are 
unphysical. The shapes of input and output functions for models satisfying this bound have qualitatively the same properties 
as shown before. From (a) we see that consistent intermediate and ultra-relativistic DBI models satisfy a very strong upper 
bound on tensors: r < 10 -15 . This limit is consistent with the bound (98) since the second plot in Fig. 11 shows that all 
consistent models with large /nl (so that (98) applies) have a blue spectrum n s > 1. 



The rest of the algorithm is identical for all priors. 

We draw the initial amplitude of the power spectrum 
of scalar density perturbations at <j) = and its value at 
the CMB scale (taken to be k C MB = 0.02 Mpc -1 ) from 
flat priors, 

P =P 3 (</) = 0) e [5 x io- 10 ,i x 10~ 8 ] 

Ps(kcMB) e [5xl0- 10 ,lxl0- 8 ]. (A3) 

We select a limited number of the slow variation pa- 
rameters in equations (82, 83) and Monte-Carlo over 
their initial conditions using the following flat priors. 





e 


[-0.05,0.1] 




Co 


e 


[-5 x 10" 3 ,5 x 10" 


-3 


3 A 


e 


[-5 x 10~ 4 ,5 x 10" 


-4 


4 A 


e 


[-5 x 10~ 5 ,5 x 10" 


-5 


5 A 


e 


[-5 x 10~ 6 ,5 x 10" 


-6 




e 


[-0.1,0] 






e 


[-0.01,0.01]. 





Other (higher order) parameters are set to zero at the 
initial conditions, and therefore remain zero for the entire 
evolution. 

With the choice of parameters in Eq. (A4), we have 
expanded H(<p) and j(<p) polynomials to 6th- and 



2nd-order in the scalar field 0, respectively. Such an 
expansion is more than adequate to capture, for exam- 
ple, the behavior of functions such as f{(j)) oc (/i + </>) -4 
through Eq. (65). We have checked that our qualitative 
results do not change by adding higher order slow 
evolution parameters to the dynamics. 



2. Evolution and Matching to CMB Scales 

We numerically evolve the flow equations forward into 
the throat (0 decreases) until the matching condition 

is satisfied. Then we need to identify the (j) value corre- 
sponding to &cmb • Note that unlike in slow roll inflation, 
the horizon exit for scalar and tensor modes for a given 
physical scale /ccmb happens at different <p values, 0^ MB 
and (j)^ MB respectively. These are related by 

a((j) s )H((j) s )>y((j) s )\cMB = a(</>t)#(0t)|cMB- (A6) 
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Rewriting this in a way to minimize numerical errors in 
the matching, 



N°-N s e (4> s )+HH(<f> s ) 7 (4> s )} 
= N° - NKfo) + In H (</*), 



(A7) 



where N® is the arbitrary number of e-folds associated 
with = 0. 

Models which never satisfy the matching condition 
(A5) are rejected. In addition, throughout the evolution, 
we impose the following constraints: 

1. In order to be compatible with the physical picture 
of single throat UV DBI inflation, the model must 
have a monotonically increasing warp factor, i.e. 
df /dcj) > must be satisfied at all times during the 
evolution. 

2. k((j)) must be a monotonic function. This requires 
e + k < 1, where 



e + k = 



2e d \n(jH) 

11 M P 



(A8) 



Since \JTpi < 1, the condition e + k < 1 implies 
dlnicsH- 1 ) 



d 



< 1 



(A9) 



M P 



In other words, the sound horizon shouldn't grow 
exponentially as the field evolves over Planckian 
distances. 

In the single throat scenarios that we are modeling, 
is constrained to be negative. To see this note 

that 



«(0) = - 



d In c s d In 7 



dN K 



dN, 



(A10) 



For a brane moving into a throat 7 = [1 — J*/) 2 ] -1 / 2 
is necessarily monotonically increasing, because 
both f((j)) and <\> are monotonically increasing as 
the brane speeds up while moving towards larger 
warping. With our convention d7V e = —Hdt this 
implies that K,(<fr) < 0. 

4. 7(0) > 1, by definition. 

5. Successful models must be able to accommodate at 
least the range of physical scales k G [1 x 10 _5 ,3] 
Mpc -1 between = and </> en d- 

Models failing any one of these requirements are rejected 
from the simulation. 



order to capture the uncertainties in the reheating en- 
ergy scale which translate into significant uncertainties 
in the number of e-folds of inflation between the CMB 
scale exiting the horizon and the end of inflation [37] 
we draw a random number of e-folds to be sampled from 
N* G [40, 70] . Then we evolve forward towards the end of 
inflation so that A^cmb — ^V en d = N*. Defining <p = e nd 
and computing / e nd = /(0end) and h end = T 3 / e nd, we 
check that h~* d G [10~ 10 , 10~ 2 ]. This is our model for in- 
flation ending via a tachyonic transition. If slow roll ends 
(i.e. e — > 1) before N* is reached, we identify this as the 
end of inflation if this occurs at least 40 e-folds from the 
CMB scale, taking the number of e-folds between 0cmb 
and e = 1 as N* . If neither of these conditions are satis- 
fied, the model is rejected. 



3. End of Inflation 



We consider two distinct scenarios for the end of infla- 
tion: tachyonic instability and the end of slow roll. In 
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